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Dirac half-metal in a triangular ferrimagnet 
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An idea is proposed for realizing a fully spin-polarized Dirac semimetal in frustrated itinerant 
magnets. We show that itinerant electrons on a triangular lattice exhibit the Dirac cone dispersion 
with half- metallic behavior in the presence of a three-sublattice ferrimagnetic order. The Dirac nodes 
have the same structure as those of graphene. By variational calculation and Monte Carlo simulation, 
we demonstrate that the ferrimagnetic order with the Dirac node spontaneously emerges in a simple 
Kondo lattice model with Ising anisotropy. The realization will be beneficial for spintronics as a 
candidate for spin-current generator. 
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Massless Dirac fermions show substantially different 
nature from ordinary electrons. The peculiar nature orig- 
inates in the characteristic energy dispersion — the nodal 
structure with linear dispersion often referred to as the 
Dirac cone. While the Dirac fermions were originally 
introduced in the relativistic quantum theory, recent dis- 
covery of graphene [HQ, a single layer sheet of graphite, 
has carved out a new direction of their study in condensed 
matter systems d, H|. In graphene, two Dirac cones ap- 
pear in the energy dispersion of tt electrons, which are 
at the K and K' points in the Brillouin zone for the 
two-dimensional honeycomb lattice. The Fermi level of 
this two dimensional conductor comes right at the nodal 
points, and the low-energy Hamiltonian is well approx- 
imated by the Weyl equation Various remarkable 
electronic and transport properties of graphene mainly 
owe to these Dirac cones in the band structure. 

The extraordinary nature of Dirac fermions in 
graphene has also attracted a great interest from appli- 
cation to electronics |3j . From the viewpoint of such po- 
tential applications, it is of great interest to control the 
characteristic band structure. Furthermore, it is also de- 
sired to control the electronic spin degree of freedom for 
the application to spintronics [6|. However, there is not 
so much flexibility in graphene, as the Dirac cone is a di- 
rect consequence of the honeycomb lattice geometry and 
the relativistic spin-orbit interaction is very weak. 

In this Letter, we propose an alternative solution for 
manipulating the spin degree of freedom by seeking pos- 
sible emergence of Dirac fermions from itinerant mag- 
nets. We show that itinerant electrons coupled to a 
well-known ferrimagnet on a triangular lattice give rise 
to the Dirac nodes in their band structure, similar to 
those of graphene. The resultant massless Dirac fermions 
are spin-polarized, and they are stable in a wide range 
of the spin-charge coupling including typical values in 
solids. We demonstrate that, by an unbiased Monte 
Carlo (MC) simulation as well as a variational calcu- 
lation, such Dirac half-metal with ferrimagnetic order 
spontaneously emerges in a minimal Kondo-lattice type 
model. The results strongly suggest the possibility of 



realizing the exotic electronic state in transition-metal 
and rare-earth compounds, which generally retain much 
higher controllable degrees of freedom than graphene. 
Such a new family will not only add a member to the 
known list of Dirac electrons in solids but also 

bring a completely new aspect by the spin polarization. 
In a half-metal, the electric current is perfectly spin- 
polarized as the low energy excitations only exist for the 
majority spin (sj. This nature works as a spin-current 
generator by filtering out the minority-spin electrons. 
Thus, our proposal opens a new frontier for the applica- 
tion of Dirac massless fermions, especially for spintron- 
ics [n|. 

Let us first discuss a naive, rather trivial approach 
to achieve a Dirac half-metal. We here consider a 
single-band ferromagnetic Kondo lattice model (double- 
exchange model) on a honeycomb or kagome lattice [see 
Figs. QJ a) andQJb)]. The model consists of the nearest- 
neighbor hopping of electrons and the exchange inter- 
action between the electron spin and localized moment, 
whose Hamiltonian is given by 

H = "* E ( C *W +H.c.) - Jj>i • Si. (1) 

Here, ^(cj^) is the annihilation (creation) operator of an 
itinerant electron with spin a =t, 4- at ith site, ai and 
represent the itinerant and localized spin, respectively; t 
is the transfer integral and J is the onsite Kondo cou- 
pling. Hereafter we take t = 1 and J > 0. 

In this model, when J is sufficiently large compared 
to the bandwidth at J = 0, a ferromagnetic order is 
stabilized by the double-exchange mechanism in a wide 
range of electron filling n = J2i a (cl a Ci a ) /2N , where N 



is the system size 
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In the ferromagnetic phase, 



the band is split into two by the large exchange coupling 
according to the spin component, and each band has ex- 
actly the same form as that for the noninteracting case 
J = 0. Hence, in principle, the Dirac half-metal arises 
for the honeycomb and kagome lattices, as the nonin- 
teracting bands on these lattices have the Dirac nodes. 
However, these situations are very difficult to realize in 



2 



(b) 



'« i< 



FIG. 1. (color online). Schematic pictures of (a) a hon- 
eycomb ferromagnet, (b) kagome ferromagnet, and (c) three- 
sublattice triangular ferrimagnet. The arrows at each site 
represent localized spins. 



solids as neither such a strong exchange interaction nor 
the honeycomb and kagome structures is easily realized 
in magnetic compounds. 

As a more realistic approach, here we propose a sim- 
ple, but rather nontrivial route to the half- metallic Dirac 
fermion systems. Let us consider the model in Eq. (pQ) on 
a triangular lattice, and the situation in which a three- 
sublattice collinear ferrimagnetic order with up-up-down 
spin configuration is realized — see Fig.QJc). By treating 
the localized moments as classical spins with |S^| = 1, 
the band structure is easily calculated by the exact diag- 
onalization of the Hamiltonian. The lower three bands 
of the totally six bands are shown in Fig. [2j the two red 
bands are of up spins, and the blue band is of down spin 
(the other upper three bands have the similar form with 
opposite spins). 

The band structure has a notable feature at the energy 
e = —J; the two up-spin bands touch with each other 
at the K and K' points in the Brillouin zone to form 
a Dirac-type point node with linear dispersion, and the 
down-spin band has the band top at the same points with 
an ordinary parabolic dispersion. See also the enlarged 
figure in Fig. Efb) and the energy dispersion along the 
symmetric lines in Fig. Efc). 

In this situation, when the electron filling is at n = 1/3, 
the two lower bands are fully occupied while the remain- 
ing bands (including the upper three) are unoccupied; the 



Fermi level is located at the nodes where the three bands 
meet. As the down-spin band has an energy gap, the half- 
metallic Dirac electrons are obtained by electron doping 
to the unoccupied up-spin band. Although hole doping 
hides the Dirac nature as the down-spin parabolic band 
is doped at the same time, the situation is avoided by in- 
troducing an additional antiferromagnetic exchange cou- 
pling between the neighboring sites, J' J2(i j) a i ' Sj Q- 
A finite J' > shifts the down-spin band to the lower 
energy and isolates the half-metallic Dirac nodes ener- 
getically, as demonstrated in Figs.[2fd) and[2^e). Hence, 
the simple ferrimagnetic order on the triangular lattice 
realizes the peculiar Dirac half- metallic state near 1/3 
filling. 

The Dirac nodes have essentially the same structures 
as those in graphene. Under the ferrimagnetic order, the 
Hamiltonian is written as 
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Here, the upper two rows correspond to the sites with the 
up localized moment (A, B sublattices) and the bottom 
row is for the down one (C sublattice) in the three-site 
unit cell. In Eq. (j2]), <j z is the z component of the Pauli 
matrix for itinerant electrons, k is the wave vector, and 
Tk is the Fourier transform of the hopping term given by 

By using 
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the k • p perturbation around the K and K r points in the 
Brillouin zone [5| and by expanding the result up to the 
first order in terms of tK x j J and tK y / J (k is the relative 
wave vector measured from K and K' points), we end up 
with the low-energy Hamiltonian which is factorized into 
two parts. One is a 2 x 2 Hamiltonian for the up-spin 
honeycomb subnetwork of the A and B sublattices, and 
the other is a localized state at the down-spin sites in the 
C sublattice. The former is given by 
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(3) 



where the sign ± corresponds to the K and K' points. 
This has an equivalent form to that of graphene. 

It is worthy to note that the Dirac nodes are formed 
immediately by switching on J. However, when J is 
small, the low-energy physics at n = 1/3 is not char- 
acterized solely by the massless Dirac fermions because 
there is a band overlap at the energy of the Dirac 
nodes. The band overlap comes from the second lower 
band for up spin, which has an energy minimum at 
k = (27r/3,0) points and its threefold symmetric points 
for small J; the minimum energy is given by e^^j = 

t/2 - y / (J + 3J / -£/2) 2 + 2£ 2 . In order for the Dirac 
nodes to be isolated at the Fermi level, this energy 
should be higher than that at the K and K' points, 
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FIG. 2. (color online). Band structures of the model in Eq. (pQ) under the three-sublattice ferrimagnetic order at J — 2. (a) 
The overall band structure of the three lower-energy bands at J' = 0, (b) the enlarged view near the Fermi level e — — J at 
n — 1/3 in the first quadrant, and (c) the cut along the symmetric lines, (d) and (e) show the results at J' — 0.05. The arrows 
indicate the spins for each band. In (a), the gray hexagon on the basal plane shows the first Brillouin zone for the magnetic 
supercell. The dashed line in (e) indicates the Fermi level in the MC simulation shown in Fig. [4l 



£k = —(J + 3J'). Hence, the Dirac nodes are energeti- 
cally isolated and play a decisive role when the condition 
(J + 3J')/t > 1 is satisfied. This condition is important 
because the necessary J and J' are much smaller than the 
noninter acting bandwidth 9£, and it is indeed satisfied in 
wide range of materials. 

So far, we assumed the presence of three-sublattice fer- 
rimagnetic order. In the following, we show that such 
order is indeed stable in the Kondo-lattice type model 
as Eq. (pQ). We here simplify the model by assuming the 
localized moments are the Ising spins taking the values 
Si = ±1. 

First, we investigate the ground state phase diagram 
near n = 1/3 by a variational calculation. We com- 
pare the ground state energy of the two-sublattice stripe 
phase and three-sublattice ferrimagnetic phase appeared 
in the previous study [15], in addition to the ferromag- 
netic phase. The results at J = 2 are shown in Fig. [3] for 
J' = and 0.05. At J' = 0, the ground state in the plot- 
ted range is dominated by the ferrimagnetic phase as well 
as the stripe phase. The different phases are separated 
by phase separation. As shown in Fig.[3jb), the introduc- 
tion of small J' largely stabilizes the ferrimagnetic phase 
near n = 1/3 as well as the stripe phase. This is because 
the itinerant electron spins are polarized parallel to the 
localized spins in the ground state, leading to an energy 
gain (loss) by the antiferromagnetic J' for the two states 
(the ferromagnetic state). 

We next examine the stability of the ferrimagnetic or- 
der at finite temperatures by an unbiased MC simulation. 
For the simulation, a standard algorithm for fermion 




FIG. 3. (color online). Ground state phase diagram obtained 
by variational calculation at (a) J' — and (b) J' — 0.05. 
The schematic picture of magnetic structure in each phase 
is shown. The white region indicates the electronic phase 
separation (PS) and the dotted vertical lines indicate n — 1/3. 
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FIG. 4. (color online). MC results for (a) the pseudo mo- 
ments M xy and \M Z \, (b) corresponding susceptibilities \xy 
and Xz-> an d ( c ) azimuth parameter ip. The data are calcu- 
lated at n = 0.34. 



systems coupled to classical fields is used [16|. In this 
method, the trace over the fermions in the partition func- 
tion is calculated by the exact diagonalization, while the 
trace over classical spin configurations is computed by 
a classical MC method using the Metropolis dynamics. 
The phase transition to ferrimagnetic phase is detected 
by using two parameter [ItJ . One is the pseudo- moment 
defined by 




(4) 



where m is the index for the three-site unit cells, and 
k) denote the three sites in the rath unit cell be- 
longing to the sublattices (A, 5,(7), respectively. We 
measure the summation M = (3/iV)J^ m S m and the 
susceptibility. The other is the azimuth parameter ip 
defined by ip = (M xy ) 3 cos 60m? where (\>m is the az- 
imuth angle of M in the xy plane and M xy = 3Mj /8 



(M 2 xy = 
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). The ferrimagnetic ordering is sig- 
2^/2/3, \M Z \ -> 1/V3, and_0 ^ 1 at 



naled by M xy 

low temperature T —> 0, respectively [HI, [Til. fl9|. 

Figure H shows the MC results at J = 2 and J' = 0.05 
in the slightly electron doped region to n = 1/3 [see also 



Fig.|2je)]. The results indicate two successive phase tran- 
sitions at T KT = 0.192(15) and at T c = 0.108(9). The 
transition temperatures are estimated by extrapolating 
the peak of susceptibilities Xxy an d Xz as N — >• oo. The 
transition at Tkt is considered as a Kosterlitz-Thouless 
type with the growth of quasi- long-range order [15[ . On 
the other hand, the phase transition at T c is a three- 
sublattice ferrimagnetic ordering. The MC result and the 
above analysis for the ground state consistently indicate 
that the three-sublattice ferrimagnetic order is stabilized 
in the vicinity of n = 1/3 in the wide range of parame- 
ters for J and J 7 , spontaneously giving rise to the Dirac 
half-metal. 

As such ferrimagnetic order was indeed observed in 
several insulating magnets [20|, [2l| , our results in the min- 
imal model will stimulate the hunt for Dirac half-metal in 
transition-metal and rare-earth compounds. The present 
results will be qualitatively robust even when extend- 
ing the model to more realistic situation. For instance, 
the ferrimagnetic state remains stable when including the 
transverse components of localized spins, at least, in the 
presence of the Ising anisotropy. Multi-band effect may 
be avoided under a particular crystal field; for instance, 
the <i-electron a± g orbital isolated by a strong trigonal 
field is a good candidate for the realization. Interlayer 
coupling, however, may open a gap at the Dirac nodes. 
Nevertheless, a straightforward stacking of layers or suf- 
ficiently isolated layers in a controlled thin film will be 
promising to preserve the massless nature. 
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